@ REPUBLiQ UE FRA NQAiSE 

INSTITUT NATIONAL 
DE LA PROPRIETE INDUSTRIELLE 



PARIS 



^ N^ de publication : 

n'utiliser que pour les 
commandes de reproduction) 

N*' d'enregistrement national 
^ int CI' : H 04 N 5/32 



2 700 909 
93 00804 



® 



DEMANDE DE BREVET D'INVENTION 



A1 



Data de depot : 27.01.93. 
(29) Priorit6 : 



uy Demandeur(s) : Sodet6 Anonyme dite: GENERAL 
^ ELECTRIC CGR — FR. 



^3) Date de la mise d disposition du public de la 
^ demande : 29.07.94 Bulletin 94/30. 

fee) Liste des documents cites dans le rapport de 
recherche prdliminaire : Se reporter^ la fin du 
present fascicule. 

feo) References ^ d'autres documents natlonaux 
apparent6s : 



lnventeur(s) : Picard Catherine, Rougto Anne et 
Trousset Yves. 



@ Titulaire(8) : 



^4) Mandataire : Cabinet Ballot-Schmit 



@) Dispositff et proc6d6 automatique de calibration g6om6trique d'un systdme d'Imagerie par rayons X. 



fen On utilise un fantdme constitu6 d'une h^lice. on mon- 
treque cette h^lice permet une calibration gtem^trique au- 
tomatique de n'importe quel syst^me d'Imagerie par rayons 
X d d^tecteur plan. 




coU^ 



I mill mil mil Hill Hill mil Hill nil nil 



DISPOSITIF ET PROCEDE AUTOMATIQUE DE CALIBRATION 
GEOMETRIQDE D'DN SYSTEMS D'IMAGERIE PAR RAYONS X 

La present e invention concerne iin dispositif et un 
precede automatique de calibration g6ometric[ue d»un 
syst^e d» imager ie a rayons X* Des systdmes de ce type 
sont utilises principalement dans le domaine medical 
5 dans le but d' analyser des structxires internes de corps 
bxmains notamment de structxires vasculaires* Le but de 
1* invention est de permettre, par une calibration aisee 
du systSme , son utilisation potar Svaluer 
quantitativement , et non plus uniquement 

10 qualitativement/ les dimensions des structures internes 
analys€es. L' invention se rapporte essentiellement aux 
systemes dont le detecteur de rayonnement est 
bidimensionnel, par exemple, sensiblement plan. Elle se 
rapporte, poxir le proced^, a\ix systemes toxirnants 

15 d'imagerie par rayons X. 

Dans le cas d'un systeme d* imager ie a rayons X 
compose d'une source de rayons X et d'un d6tecte\ir 
bidimensionnel, toumant ou non, 1* operation geom4trique 
intervenant dems la production de 1» image est une 

20 projection conique d»\in objet analyse, deploy^ dzms un 
espace S trois dimensions (3D) , sur un espace §l detix 
dimensions (2D) qui est celui du plan de projection 
correspondant au plan de detection. Les pcirametres 
geometriques dScrivant complStement la projection 

25 conique sont au nombre de neuf . Ces neuf paramStres 
peuvent etre pr^sentes selon diffSrentes 

parametrisations . 

Une des parametrisations possibles est la suivante: 
- trois premiers parametres concernent les trois 

30 coordonnees de la source S de rayons X dans le repfere 



R de I'espace 3D dans lequel est place I'objet 
analyst; 

_ trois autres pEurametres concernent les trois angles 
d'Euler du changement de repSre, entre le repere R de 
I'espace 3D de I'objet et un repere R' liS ^ la 
projection conique, c'est a dire formS d'un rep&re 2D. 
du plan de projection P et d'un troisieme axe 
orthogonal au plan de projection; 

- trois derniers parametres concernent les trois 
coordonnees d'un point arbitraire I du plan de 
projection P (par exemple le centre de 1' image) dans 
le repere R. 

Pour simplifier les calculs dans 1* invention on 
verra par la suite qu'il y a lieu de preference de 
retenir des coordonnfees cylindrigues pour les trois 
premieres et trois demieres coordonnees. Le repSre R 
est alors un repere cylindrique- 

La connaissance de tout ou partie de ces paramStres 
est tr&s souvent utile en radiologie, notamment 
lorsqu'une donnSe quantitative sur un objet 3D est 
estimSe a partir d'une mesure faite dans la projection 
2D, ou dans plusieurs projections 2D, de cet objet. Des 
paramStres de chaque point de vue doivent alors' etre 
connus. Un point de vue conceme 1 » orientation du 
systSme par rapport 4 I'objet. A titre d' exemple, on 
retiendra le problSme classique de 1' estimation de la 
section d'un vaisseau aj partir d'une ou plusieurs 
projections 2D. Or il est souvent impossible, ou trop 
imprecis, d'accSder a ces parametres de fa^on directe, 
c'est a dire, par exemple, en mesoirant directement sur 
le systeme d' acquisition la distance entre la source de 
rayons X et le dStecteur. 

On appelle calibration g^omStrique d'un systdme 
d' imager ie, 1' operation qui aboutit a la connaissance 
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indirecte precise des parametres geometriques qui 
interviennent dans la production d'xane image. Le 
principe, classique, est base sur 1* utilisation d*\in 
fantome geometrique connu dans I'espace 3D, et dont on 
acquiert la projection 2D. La suite des operations 
effectuees dans ce but comporte les etapes suivantes: 

- on dispose d'un objet connu, un fantdme de 
calibration prSsentant un certain nombre de points 
caracteristiques dont la position dans l*espace est 
connue par des coordonnees mesurees par rapport a xm 
repere propre d cet objet; 

- on acquiert 1» image de ce fantome dans les 
conditions gSomStriqoies d*un point de vue (ou 
incidence) que l»on veut calibrer; 

- on reconnalt les projections des points 
caracteristiques dans 1« image. Pour cela on associe 
ctiaque point caracteristique de 1» objet a sa trace 
dans 1' image acquise, la projection; 

- on inverse au sens mathematique le systeme 
d* equations decrivant la projection; 

- et on obtient finalement 1» ensemble des paramStres 
de la projection pour le point de vue donne. 

Dans I'etat de I'cirt, la reconnaissance des points 
caracteristiques dans la projection du fantSme de 
calibration, est faite "manuellement" par un opSrateur 
bumain. Par exemple, une forme de fantome de calibration 
geometrique souvent utilisee est celle d'un cube, a\ix 
buit coins duquel sont disposees des billes mStalliques 
opaques aux rayons X. Parf ois des billes supplementaires 
sont rajoutees pour augment er la precision geometrique 
de la calibration. Or, la projection au sens des rayons 
X produit une image "en trsuisparence" de 1» objet, dans 
laquelle, selon 1 • orientation du cube, il peut etre tr&s 
difficile d'associer sans se tromper les traces 2D des 



4 

billes avec lexirs homologues 3D. On peut se refferer, 
pour I'etat de I'axrt de la calibration geometrique en 
radiologie axix articles suivants: 

(1) D.L. Parker, J^Wu, D.L Pope, Van Bree, G.R. 
5 Caputo and H.W. Marshall, "Three-dimensional 

reconstruction and flow measurements of coroncury 
arteries using multi view digital angiography", in New 
Developments in Quantitative CoronsLry Arteriography, 
J.C. Reiber and P.W. Serruys Eds., pp. 225-247, Kluwer 
10 Academic Publishers, 1988; 

(2) D.J. Hawkes, A.C.F. Colchester and C.R. Mol; "The 
acciirate 3-D reconstruction of the geometric 
configuration of the vascular trees from X-ray 
recordings," in Physics and Engineering of Medical 

15 Imaging, R. Guzzardi Ed., Nijhoff, 1987. 

(3) M. Garreau, J7I* Coatrieux, R. Col lor ec and ^C. 
Chardenon, "A knowledge-based approach for 3-D 
reconstruction and labeling of vascular networks from 
biplane angiographic projections", IEEE Medical 

20 Imaging, vol. 10, N*^2, pp. 122-131, June 1991. 

Cette necessite d » intervention humaine est un 
inconvenient majeur, et peut mSme devenir prohibitive 
dans certains cas, lorsque le nombre de points de vue Sl 
calibrer est importcuit. C*est le cas en particulier pour 

25 tout systeme d* acquisition de projections 

bidimensionnelles par rayons X dans le but d'une 
reconstruction 3D. Un tel syst^e de reconstruction est 
par exemple decrit dans la demande de brevet frangais 
N*» 2 644 590 deposSe le 20 mars 1989. Quelle que soit la 

30 m€thode de reconstruction 3D, il est necessaire au 
prSalable de connaltre parfaitement les paramStres 
geometriques qui caracterisent chacune des projections. 
Dans le cas probable oil une m^thode de calibration 
indirecte est necessaire, il semble indispensable 



d'effectuer cette calibration de fagon automatique. 

La presents invention propose done un dispositif 
permettant xine calibration automatique, qui s' applique Sl 
tout systeme d^imagerie bidimensionnelle a rayons X, et 
5 en particulier Si tout syst&ae utilise S des fins de 
reconstruction 3D. 

Le dispositif de 1* invention est tel que la 
reconnaissance, sans equivoque, des correspondances 
entre les billes et leurs traces dans les projections 

10 est automatique* Ceci est obtenu en choisissant un 
fantdme dans lequel les billes sont distributes, de 
procbe en proche, en une succession telle que des 
altitudes de billes, mesurees le long de I'axe de 
rotation du systeme d* imager ie, et siirtout un axe du 

15 fant5me, soient monotones croissantes (ou decroissantes) 
avec un numfero d'ordre des billes dans la succession. 

Dans ce but 1» invention a done potir objet xin 
systeme d» imager ie par rayons X avec xin dfitectettr 
bidimensionnel situe en vis a vis d"\in tube a rayons X, 

20 ce dttecteur et ce tube pouvant dans certains cas 
tourner ensemble autour d*un axe, ce systSme d» imager ie 
par rayons X comportant un dispositif de calibration 
geometrique muni d'xin fantSme de dimensions connues et 
de moyens de mesure pour deduire, S partir d* images en 

25 projection du feuitame sur le detecteur bidimensionnel, 
des coefficients de calibration de ce systSme d'imagerie 
par rayons X, caracterise en ce que le fantSme comporte 
une succession de proche en proche de structures 
cellulaires a eibsorption radiologique difftrente, plus 

30 elevee ou plus faible, par rapport a levir environnement, 
chaque structure cellulaire etant identifiable 
automatiquement par une caractSristique ordonnte 
d' identification de cette structure cellulaire. 

Un autre probleme a resoudre poxor 1' invention est 



lie au nombre des calibrations S entrepr endr e . Un 
systeme <!• imager ie toxirnant de ce type est en effet 
susceptible d'occuper 1000 positions angulaires, ou 
plus, sur un tour* Des calibrations pour toutes les 
5 positions angulaires, meme men€es automat iquement, sont^ 
trop longues a ef f ectuer povir pouvoir Stre r^petees 
periodiquement, toutes les semaines ou tous les mois. II 
faut en fait prSvoir une journ^e de travail. 

Le proc^dg de 1» invention vise S rSsoudre cet autre 

10 problSme en faisant 1 » bypothSse que, si les dimensions 
tbSoriques de construction du systSme d ' imagerie ne 
peuveiit pas etre considSrees comme respectees dsuis le 
systeme rSel, on peut neemmoins admettre que le systdme 
d' imagerie est Sl rotation circulaire parfaite ou presque 

15 parfaite. L'idSe du procedfe est alors de recbercher- les 
parametres intrinseques du systeme: en definitive cexix 
qui donnent en plus la direction de son coce de rotation. 
A partir de ces parametres intr insecfues , cju^on peut 
connaltre trSs prScisement, on peut ensuite calculer 

20 pour cbaque position angulaire des coefficients de 
calibration. On montre que, par exemple, xine centaine de 
calibrations conduisent a une connaissance tres precise 
de ces parametres intrinseques. Avec cette connaissance, 
on calcule 900 autres "calibrations". On en connalt 

25 alors 1000 au total en n"en ayant mesurees que 100. Le 
gain de temps est considerable. 

Ii» invention a done egalement pour objet un precede 
de calibration geometrique d'xin systdme d' imagerie pear 
rayons X H dStecteiir bidimensionnel tournant autour d»un 

30 axe, en vis a vis d'un tube a rayon X et comport cuit les 
etapes suivantes: 

- on place un fantome de dimensions connues entre le 
t\abe et le detecteur, 

- on mesure, pour une position en rotation de 



1' ensemble txibe detecteiir par rapport au fcuit6me, dans 
1 * image radiologique du f antome pro jetee sur le 
detectevir bldlmensionnel, des coordonnees de lleux 
images de points caracteristiques du f antome , 
5 - on en deduit, poiir cette position en rotation de 

1' ensemble tube d^tecteur, et dans un repere associe k 
la position du f antome, des coefficients de calibration 
relatifs aux positions respectives d'un foyer de 
rayonnement du tube Si rayons X et du detecteur 
10 bidimensionnel , 

- et on reit&re ces deux dernieres etapes pour des 
positions desirees diffSrentes de 1' ensemble tube 
detecteur r 

caracteris^ en ce que 

15 - on traite les coefficients de calibration relatifs a 

un certain nombre de positions en rotation de 
1» ensemble tube detecteur pour en extraire des 
parametres intrinseques du systeme d'imagerie psir 
rayons X qui sont indSpendants de la position en 

20 rotation de 1» ensemble tube detecteur, 

- et on deduit de ces parametres intrinseques des 
coefficients de calibration relatifs H une position 
quelconque en rotation de 1' ensemble tube dStecteur en 
evaluant les valeurs de f onctions analytiques de ces 

25 param&tres intrins&ques pour une valeur angulaire de 

position en rotation de 1» ensemble tube d€tecteur. 

invention sera mieux comprise a la lecture de la 
description qui suit et & I'examen des figxires qui 
I'accompagnent. Celles-ci ne sont donnees qu^S titre 
30 indicatif et nullement limitatif de 1* invention. Les 
figures montrent: 

- Figure 1: un systeme d» imager ie utilisable pour 
mettre en oeuvre 1 • invention avec la presentation des 
differents reperes d» acquisition; 



- Figiiires 2a S 2d: respectivement des vues de prof il 
et de face du fcintome selon 1» invention/ la trace des 
bilLes dans les projections, et un detail de mise en 
place des billes peur inrarustation; 

5 - Figxire 3: la mise en place du famtSme dans le 

systeme d» imager ie; 

- Figxires 4a et 4b: des images en projection d'une 
forme particuliere du fantdme apres traitement; 

- Figures 5a a 5d: des representations geometriques 
10 permettant tine meilleure comprehension du proced4; 

La figure 1 montre \in systeme d'imagerie utilisable 
pour mettre en oeuvre 1* invent ion • Ce syst&ae comporte 
un tube 1 S rayons X avec un foyer S de rayonnement, Un 
detecteur 2 de rayonnement , par exemple du type 

15 amplif icatexir de lumineuice utilise en radiologie, est 
place en vis a vis du tube 1» Une normale au plan de 
detection du detecteur 2 passe sensiblement par- le 
foyer S et le centre I de ce detecteur. Ceci peut atre 
une necessite pour la formation d' images classiques, 

20 mais n'est pas directement lie au problSme de la 
calibration. En psirticulier, on peut calibrer, avec 
1' invention r vin tomographe du type a translation 
liomotlaetic[ue. L» ensemble peut toumer avec xon mecanisme 
(non represents) de type connu pour les 

25 tomodensitomfetres, ou avec un arceau, autoxir d»un axe 3 
montr€ ici verticalement mais qui, dans la pratique, est 
normalement horizontal. Dans les cas des tomographes, il 
n'y a que des translations et un centre d'tiomothStie, 
sans mSme nScessairement un axe. II en est de meme bien 

30 sflr des appareils staticfues. Le systdme comporte 
egalement un processeur de traitement 4 qui fait subir 
aux signaux Slectriques mesxxres par le detecteur 2 im 
traitement contenu dans un programme 5 stocke dans une 
memo ire programme 6. Des images, ou des resultats 



statistiques resulteuit de ces traitements, peuvent etre 
visualises sur xin monitexir ?• 

Ij ' invention permet de rSsoudre le probleme de la 
reconnaissance automatique des points Iiomologues 2D/3D, 
5 en utilisant comme fantome de calibration im fajitome 
particulierement adapts montre sur les figures 2a a 2d 
et 3. 

Ce fantome est dans un exemple, un cylindre 8 de 
plexiglass, qui sert de support a tme s&rie de petites 

10 billes metalliques incrustees dans la parol du cylindre, 
et disposees rSgulidrement selon une helice. Les billes 
ont toutes un diamStre identique, par exeiaple 1 mm, sauf 
la bille du milieu, qui a un diamStre plus Sieve que les 
autres, par exemple 3 mm* Les paramStres de I'hfelice 

15 sont grossierement adaptes au systSme d* acquisition des 
images, de fa9on a ce que sa projection soit la plus 
grande possible sur le detecteur tout en restant 
contenue dans 1' image. De m€me, le pas de I'hSlice et le 
pas d'echcintillonnage des billes sur I'helice sont 

20 clioisis de telle sorte qu'un nombre suffisant de billes 
apparaisse dans la projection, sans qu'elles puissent se 
superposer. Les figures la Id sont cotees en 

millimetres. Pour eviter toute deformation du fantome 8, 
notamment par vieillissement du plexiglass le cylindre 8 

25 est obturS des deux cotes pax un couvercle 9 qui en 
assvire la rigidity. 

Le fant6me 8 est adapte par exemple a ^xn systSme 
d' acquisition constituS d'une soiirce de rayons X et d'un 
detecteur de champ rond de diametre 37 cm, situS a une 

30 distance d» environ 135 cm, alors que I'objet est a 
environ 100 cm du foyer S. Ces chiffres sont des 
exemples qui ne sont en rien restrictifs, de meme que 
les dimensions du fantSme indiquees en annexe. Dans 
1' exemple il y a 61 billes. Une des coordonnees des 



billes, 1* altitude mesuree sur xin axe longitudinal Z, 
evolue de proche en proche de 5 mm en 5 mm avec une 
precision de 5/lOOeme. II y a 20 billes rSparties 
regulierement par tour de I'helice. 
5 Les trois particulctrites principales du fant6me 8 

sont les suivantes^ D'abord le fantdme 8 dispose d'un 
nombre suffiscoit de billes pour la calibration, Elles 
sont bien r^parties dans I'espace 3D. Leurs projections 
sont egalement bien reparties dans le plan P (condition 
10 sine qua non pour une calibration 3D) . La rSgularite de 
la distribution n*est cependant pas indispensable. Dans 
ce but les billes peuvent se nxmeroter facilement en 
fonction de leur cote selon I'axe Z. Ceci est obtenu ici 
selon une formule trSs simple: 
15 2 = h.n+zO, 

z = coordonnees selon Z, 
zO = coordonnee de la bille numero 0, 
h = pas d'ecbantillonnage des billes selon Z, 
n = niomero de la bille. 
20 C'est en fait la relation d'ordre qui nous importe ici: 

si zl>z2^ alors nl>n2 
zl^ z2 = coordonnees des billes 1 et 2 selon Z 

nl, n2 = numero des billes. 
Si le fantSme est positionne de fagon judicieuse 
25 devant le d^tecteur de rayons X, c'est a dire 
essentiellement de telle fa^on que I'axe Z soit a peu 
prSs parallele au plan P du detecteur 2, et si la 
conicite du faisceau n"est pas trop importante, cette 
relation d'ordre se retrouve egalement dans les 
30 projections 2-D. 

En effet, figure 3, soit V l»axe dans le plcui P de 
projection qui correspond a la projection de I'axe Z. On 
aura la relation suivante: 

si vl>v2, alors nl>n2 



vl, v2 = coordonnees des projections 
des billes 1 et2 selon V 
nl, n2 = nximeros des billes 
Dans la suite ^ et par souci de simplicite, on 
5 supposera c[ue I'axe V coincide sensiblement avec I'cixe 
vertical dans les images de projection, c'est & dire est 
parallele aux colonnes de 1* image. Si ce n'est pas le 
cas, on peut introduire une rotation dans le plan de 
projection po\ir s»y ramener. 

10 Toujotirs par souci de simplicite poiar la suite, on 

definit le repSre dans l*espace 3D comme etant lie S 
I'helice: c'est h dire que I'axe Z de l»helice d6finit 
tin des axes du repere R, les deux autres axes Stcmt 
perpendiculaires 3. Z. Dans ce repere, les coordonnees 

15 des billes s'expriment tres facilement. Ceci signifie 
gue,^ une fois la calibration faite, le systSme possSde 
\in axe vertical qui est celui qu'occupait I'axe du 
fantdme au cours de la calibration. Des reconstructions 
3D relatives S. des calibrations diff^rentes ne sont done 

20 pas directement comparables. Ceci sera par centre 
possible avec le precede de 1» invention vu plus loin. On 
a represente dans la figure 3 le fantome 8 positionne 
entre la soxirce S de rayons X, et le detecteur D de 
rayons X. Le fantome est positionn^ de telle fagon que 

25 son cixe principal Z est a peu pr^s parallele au plan 

et tel q[u"il se projette selon I'axe V du detecteur. . 

Enf in les billes peuvent se nxmeroter de fagon 
absolue. K&sie si le haut et le bas de la projection de 
l»h61ice sont tronquees (c'est a dire que l»on ne voit 

30 pas la premiere bille ni la demi&re bille du fantome) , 
on a xine reference absolue c[ui est la grosse bille 
particuliere du milieu. La bille particuliere, ou deux 
billes particulidres , ou autres, sont done de preference 
placees dans le fantome a un endroit qui a de grandes 
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chances d'etre dans le champ: au milieu du fantome. 

Poxir qu'on comprenne le principe mis en oeuvre dans 
le precede de calibration, on a decrit dans les 
paragraphes precedent un cas particulier de f cmt6me qui 
repond aux criteres exposes. En fait, tout fantSme ayant 
les particularites suivantes pearmet une calibration 
automatique, Ces pcirticularites sont: 

- Le fant6me doit disposer de structures cellulaires, 
pour simplifier des billes, visibles en rayons X, et 
bien reparties dans I'espace 3D. Leur contraste 
radiologique par rappoirt a leur environnement doit 
Stre eleve. Les projections des billes doivent Stre 
contenues dans 1» image 2D au mo ins pour les billes 
centrales. II n'est pas du tout nScessaires que, le 
support des billes soit un cylindre creux z de 
plexiglass. 

- II doit exister une relation d'ordre entre : un 
numero affectable aiix billes et leurs coordonnees, 
aussi bien dans I'espace 3D selon I'axe Z, que dans la 
projection selon tui axe, peir exemple I'axe V, . et 
quelle que soit 1 ' or ientation du fant6me dans 
I'appareil. L'axe du fantome doit Stre sensiblement 
parallele au plan du detecteur. La solution des billes 
reparties reguliSrement sur une helice comme dans le 
cas particulier du fantome decrit, n'est pas du tout 
la seule. 

- On doit pouvoir reconnaltre au moins xine des billes 
parmi les autres, de fa9on Sl pemnettre une 
nximerotation absolue meme si on ne voit pas la 
premiere et la demiere bille dans la projection. 
Cette bille de r^f^rence peut se distinguer par autre 
chose qpie par son diam^tre, par exemple sa forme, son 
opacite aux rayons X, son absence, ou encore xine 
irr^gularite du pas de distribution, ou autre. 
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A titre d" exemple,- on citera line autre realisation 
possible pour le fantome: les billes sont reparties sur 
une hSlice "degeneree" , qui est portee par une sphere de 
plexiglass au lieu d'Stre portee par un cylindre. Dans 
5 toute la suite, on raisonne avec le fantome pr§fer6 
"Ixelice" tel que decrit ci-dessus. Toute la mfethode 
exposSe ci-apres s» applique egalement avec un autre 
fajitdme rSpondant aux criteres. 

On indique ci-apres la suite des operations a 
10 effectuer pour calibrer le systeme. 

1. Le fantome 8 decrit ci-dessus est positionnS de 
telle faqjon que son axe principal (axe Z) soit a peu 
pres parall&le au plan du detecteur. C'est le cas pcur 
exemple lorsqu'on le positionne couche sxir une table 
15 d'examen d»\in syst&me angiographique, et que I'arceau 

qui porte le tube a rayons X et le detecteur, est 
positionne dans \in plan perpendiculaire a la table. 
Dans cette position, les projections des billes 
verifient la relation d'ordre ci-dessus. 
20 2. Une projection numerisee du fantome est acquise 

dans les conditions d» incidence que I'on veut 
calibrer. 

3. L» image est transferee au processeur 4 de 
traitement d» image. 
25 4. L« image est corrigee des distorsions geometriques 

eventuelles, par exemple avec une methode telle que 
celle decrite dans la demande de brevet frangais 
N** 2 633 793 deposSe le ler juillet 1988. 

5. L» image corrigSe I est ensuite soiimise a un 
30 algorithme de detection, localisation, et 

labellisation automatique des billes decrit ci-aprSs. 

6. Les coordonnSes des billes dans le repere 3D, R, et 
de leurs projections dans le repere 2D du plan P, sont 
exploitees pour estimer les parametres geometriques de 
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la projection: les coefficients de calibration. Un 
algorithme base s\ir une minimisation d*erreur par ^la 
technique du gradient conjugue est propose et 
egalement dScrit ci-apres dans ce but. 
5 On decrit ici un exemple possible de l"algorithme 

de calibration^ qui permet de passer de 1* image en 
projection du fcintome de calibration, figure 
4a - image I, a des parametres geometriques intrinsdques 
de la machine. On montre tout au long de cette 
10 description les images resultant des diffSrents 
traitements . 

1* D^abord on segmente 1» image I. L'objectif est ici 
de supprimer dans 1» image I le signal de la projection 
du cylindre de plexiglass qui supporte les billes. On 

15 utilise pour cela des outils de morphologie 

mathematique: transformation chapeau haut de forme, 
avec un Element structurant de forme ronde et „de 
taille superiexire a la taille des billes dans 
1» image 1. On obtient une image SI qui ne contient 

20 plus que le signal correspondant a 1' absorption des 

billes. Les outils de morphologie mathematique sont.de 
type connu et publies dans de nombreux ouvrages. 

2. L» image SI doit etre purifiee des artefacts 
residuels, provenant des bords du cylindre. On utilise 

25 pour cela de nouveau une transformation de morphologie 

mathematique viseuat a supprimer les structures 
verticales dans 1» image. L* image r^sultat est une 
image S2 visible figure 4b en negatif. L' image S2 
contient done toutes les billes (signal positif ) sur 

30 un fond a valeur nulle, au bruit pres. 

3. On segmente ensuite 1» image S2. L'objectif est ici 
de ne garder que le signal provenant de la grosse 
bille centrale du fantome. On utilise po\ir cela xine 
erosion de 1» image S2, d'une taille superiexire a la 
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taille des petites billes mais infSriexire a la taille 
de la grosse bille. On obtient une image S3, qui ne 
contient que la grosse bille. 

4. On detecte ensuit:e toutes les billes dans 1* image 
S2. On veut ici localiser au pixel pres les 
projections de toutes les billes dans 1' image. On 
commence Sl la ligne 0 de 1» image (en haut du 
detecteur ) . On parcourt toute 1 • image du haut vers le 
bas. Un signal positif rencontre dans I'dlmage est 
dSclarS etre lane bille si ce signal est superieur Sl tm 
certain seuil 1. On donne le numero 0 a la premiSre 
bille rencontree,. puis le numSro est incremente a 
cbaque bille rencontree^ jusqu'en bas de 1' image. On 
est stir que I'ordre des numeros est le bon, grSce aux 
proprietes de 1'h.elice decrite ci-dessus. 

5. On calcule ensuite les coordonnSes exactes des 
billes. I**objectif ici est de connaltre les 
coordonn^es de chacune des billes dans 1» image, et 
ceci avec une precision siab-pixel. Pour chacpie bille 
dans 1* image S2, on effectue un calcul de centre de 
gravite en niveaux de gris de 1» ensemble des pixels 
autour de la bille d^tectee, et tels que leur niveau 
de gris soit superieur a un deuxieme seuil, seuil T2 
(seuil T2 < seuil Tl) . Un calcul de ce type est par 
exemple entrepris dams la deuxidme demande de brevet 
frangais citSe ci-dessus. Son principe repose sur le 
fait que les billes sont suffisamment fines pour ne 
pas absorber tout le rayonnement X. Alors, en 
1* absence de saturation d' absorption line precision 
infer ieure au pixel d* image est obtenue. 

6. On detecte ensuite la grosse bille dans 1» image S3; 
on connalt alors ses coordonnees au pixel pres. On 
compare ses coordonnees avec celles de toutes les 
billes d&tectSes dcins 1» image S2, et on reconnalt 



16 

alnsi dans la liste de toutes le bille celle qui est 
la grosse bille. 

7. La grosse bille etant reconnue dans la projection, 
on peut maintenant renumeroter toutes les billes de 
5 fagon absolue, et ainsi etre assure d'une 

correspondance exacte, pour ctiac[ue bille, entre son 
nxamero et ses coordonnees 3D, connues par ailleurs, et 
les coordonnSes 2D de sa projection. 

Connaissant, pour chaque bille, ses coordonnees 3D 
10 et les coordonnees 2D de sa projection, on peut calculer 
l^errexir de reprojection associ€e Sl un choix de valeiirs 
arbitraires pour les 9 parametres de la projection 
conique. Cette erreur est obtenue par la difference 
entre la position calculSe avec ces valeurs arbitraires 
15 et la position trouvee dans 1» image du fant6me, de 
calibration. On peut done estimer les paramStres, ou 
coefficients de calibration associSs au point de vue,. a 
la position en rotation, que l»on veut calibrer en 
choisissant les parametres minimisant un critere 
20 d» erreur de reprojection pour 1» ensemble des billes 
detectSes. 

L'algorithme propose d' estimation des parametres 
utilise la technique de minimisation du gradient 
conjugue appliquS au critSre d' erreur quadrat ique 

25 moyenne. Cette technique class ique necessite le calcul 
du critSre et de ses dferivees premieres en fonction des 
paramStres. Cet algorithme est publie aux noms de W.H. 
Press, B.P. Flannery, S.A. Texikolsky cind W.T. 
Vetterling," Numerical Recipes in C", Cambridge 

30 University Press, 1988, chapitre 10.6, page 317 et 
suivcmtes. D'autres algorithmes, decrits au meme endroit 
dans cet ouvrage, et denommes Simplex et Newton sont 
aussi utilisables. Le premier. Simplex, t l*avantage de 
s* adapter Sl n'importe quel systSme coherent de 
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parametrisation. II ne converge cependant pas 
rapidement. L • algor ithme du gradient conjugue converge 
plus vl-te mals necessite xm choix de parauiiet:risa1:lon t:el 
qoie les expressions analytiques des paramStres soient: 
5 derivables. CI est le cas de la paramfetrisation pr6sent6e 
ci-dessus. Les calculs relatifs a cet algorithme on ete 
prSsentes au congrds "3D Advanced Imaging Processing in 
Kedicine" d Rennes en France en 1992 au coxirs de la 
14eme conference cuinuelle IEEE EMBS, du 2 au 4 Novembre 
10 1992, par Anne Rougee, Catherine Picard, Cyril Ponchut 
et Yves Trousset. 

Les etapes en sont les suivantes: 
1. Initialisation: les paramdtres doivent §tre 
initialises a des valeurs realistes, ce qui permet de 
15 lancer la minimisation dans une region proche de la 

solution. 

2* Normalisation: autour de cette valeur initiale on 
peut calculer la valeur d'un critSre cozrrespondant S. 
un deplacement d'xine unitS sur chaque axe de l*espace 

20 des paramStres, puis effectuer un diemgement de 

variable dans cet espace par changement d'echelle sur 
chague cixe af in que tout deplacement unitaire stir 
chaque axe provoque le meme effet sur le critere. 
3. Minimisation: partant de la valeur initiale, et 

25 tenant compte d'un changement d'echelle trouve au 

coxirs de l^etape de normalisation, on applique la 
technique iterative du gradient conjuguS poxir 
minimiser l*erreiir quadratique moyenne sur 1* ensemble 
des billes detectSes. A chaque iteration on calcule 

30 une nouvelle estimation des paramdtres en fonction des 

derives du critSre. Le processus est stoppg lorsque 
I'ecart relatif entre 1* estimation courante et la 
suivante est infer iexire a un seuil donnS. 

On appelle machine de reconstruction 3D une machine 
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qui effectue des radiographies d'\m objet sous 
differentes incidences obtenues en tournant autoxir d'un 
axe. Ces radiographies sont ensuite utilisees pour 
reconstruire une information 3D. Cette machine peut pcir 
5 exemple etre construite S. base d'un anneau rigide sur 
lequel tourne !• ensemble tube detecteur auquel cas l"axe 
de rotation est fixe et ne change jamais au cours du 
temps. Ou bien elle peut etre construite a base d'un 
arceau inclinable supportant 1* ensemble tube dStecteoir^ 

10 auquel cas l"axe de rotation peut etre positionne de 
diffSrentes fagons dans I'espace. 

Plutot que de calculer potir chaque point de vue, 
pour chaque position en rotation de 1' ensemble tube 
detecteur, les param^tres dits coef f icients- de 

15 calibration, on peut selon un proced6 de 1» invention, 
calculer des paramStres intrinseques du^ syst&ae 
tournant. Ces paramdtres intrinsSques sont alors des 
fonctions smalytiques de la position en rotation de 
1« ensemble toxirnant. II suffit ensuite, pour obtenir les 

20 parametres de calibration relatifs a une position de 
calculer la valeur des paramfetres intrinsSques pour une 
valetar donnee d • or ientation du systSme tournant. Les 
systdmes de I'art connu delivrent d4ja normalement une 
mesure, le plus souvent sous forme d'un signal 

25 electrique, representcint leur orientation dcins l*espace. 
II suffit d'utiliser cette mesure pour obtenir les 
coefficients de calibration relatifs a cette orientation 
en la portant dans les expressions analytiques des 
paramStres intrinsdques . 

30 Poxir estimer les param&tres intrinseques dScrivant 

sous forme condensee la gSometrie globale d'tane machine 
de reconstruction 3D, on peut proceder de la fagon 
suivante: estimer les parametres de la projection 
conique, vue par vue, lors d'une premiere etape, puis 
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estimear les pareunetres intrinseques a partir de 
1* ensemble de ces parame-tres lors d'xxne dexixieme etrape. 
Piremidre etape 

Dans le cas d'une machine de reconstruction 3D^ la 
5 methode de calibration automatique est particuliSrement 
bien adaptee pour estimer les coefficients de 
calibration vue pear vue au cours d"iine rotation. 

On positioime I'helice siir la table d'examen, son 
axe Z approximativement centre svir I'axe de rotation de 
10 la machine. Pour chacun des points de vue d' acquisition 
de la machine au cours d'une rotation, on va acquSrir la 
projection de I'helice, et traiter ainsi toutes les 
images acquises, avec la methode d^crite prec^demment. 

La parametrisation proposee - avec des coordonn6es 
15 3D cylindriques - est adaptSe a 1» estimation des 
paramStres le long d'une trajectoire circulaire autour 
de 1 • axe Z . Les etapes d ' initialisation et de 
normalisation sont effectuSes pour la premiSre vue 
seulement. Les valeurs des pciramdtres trouves pour la 
20 premiere vue sont ensuite utilis€es pour initialiser 
1» estimation des pareuaStres de la vue suivante et ainsi 
de suite jusqu'^ la dernier e vue. Ceci permet 
d'accelSrer les calculs. 

On obtient done les paramdtres geometriques povir 
25 chaque point de vue, qui sont stockes stir le disque 
num&rique de la machine, et serviront dans toute gtape 
de reconstruction 3D. 
Deuxidme &tape 

II est possible selon 1' invention de passer de la 
30 connaissance des coefficients de calibration, vue paor 
vue, de la machine, a la connaissance de pareimetres 
intrinseques dScrivant sous forme condensee la gSomfetrie 
globale de la machine. 

Soit N le nombre de points de vues calibres. Par 



exemple N vaut 100. Dans le cas oCl chaque point de vue 
est calibre independamment de son voisin, on dispose de 
9N parametres poxir decrire la geometrie des 
acquisitions • Or dans le cas d'un systeme tube dStecteur 
5 lies, qui effectue une rotation autotir d'un axe fixe, on 
voit bien que ces 9N parametres sont redondants. La 
geometrie d • acquisition de la machine de reconstruction 
3D peut se parametrer avec iin nombre beaucoup plus 
limits de parcaaStres , qui sont; de preference: 
10 - les coordonnees d»un isocentre O du systSme 

(3 paramdtres) ; 

- la direction de I'axe Z» de rotation passant par 
1* isocentre O. Elle est donnee par 2 angles 
(2 parametres) ; 

15 - les rayons rj et rs des trajectoires de I et de la 

soxirce S (2 pareimetres) ; 

- La difference de cote (suivant I'axe de rotation 
Z') entre la trajectoire de la source et celle de I 
(1 paramStre) ; 

20 - I'ecart angulaire e entre 1' angle polaire de S par 

rapport a Z' et 1' angle polaire de I par rapport a Z' 
(1 paramStre) ; 

- les 2 tilts cci et a2 du detectetir, les angles 
d • orientation du dStecteur par rapport a l*axe Z» 

25 (2 param&tres) ; 

- l^amgle a3 entre les colonnes de 1» image et la 
projection de l"axe Z* dans 1» image (1 parametre) 

On obtient done un total de 12 paramdtres. D'autres 
choix pour ces douze param&tres sont possibles. Celui-ci 
30 est prSfSrS pour des considerations analogues aux 
precedentes: simplification des calculs. On notera que 
dans ce cas le referent iel dans lequel sont exprimees 
les valeurs de ces parcimetres intrinseques est 
independant de la position du fajitome. II est done 



2700909 



21 

transposable d'une calibration a 1' autre. Une methode 
pour estimer ces 12 parametres en fonction des 9N 
pcarametres de depart peut atre par exemple la suivante: 
1 • Trouver la direction de 1 ' axe de rotation 2 » : 

On considdre les N points In (n=l,N) formant la 
trajectoire de I. En theorie, cette trajectoire est 
plane et ortliogonale a Z • . On ctierche done le plan le 
plus proctie, au sens des moindres carrSs, de 1> ensemble 
des In on en d^duit la direction de Z • , qui est la 
normale ^ ce plan. 

Soient (Xn,Ynr Zn) les coordonnees d»un point In- 
On cherche le plan d* equation z = ax + by + c qui soit 
le "plus proctie" de 1' ensemble des In* On appelle 
d(a,b,c) la distance entre le plan z = ax + by + c et 
1* ensemble des points In. On rectierclie les coefficients 
(a,b,c) qui minimisent la distance d(a,b,c) . On 
d6f init: 

dCa,b,c) = S (axn + byn + c - Zn) ^ 
et on trouve alors tr%s facilement le triplet (a,b,c) 
qui minimise d(a^b,c) par regression lin^aire. 

2. Trouver les coordonnees de l»isocentre, et le rayon 
de la trajectoire de I: 

L'isocentre est le centre a la fois de la 
trajectoire de I et de celle de S. On va done rechercher 
ce centre en analysaoit la trajectoire de I. 

Remarque: dans cette etape, on cherche les 
parametres (centre et rayon) de la trajectoire de I, 
puis dans l»€tape suivante on cherchera le rayon de la 
trajectoire de S. On aurait pu, en toute logique, faire 
1' inverse, c^est 3. dire commencer par chercher le centre 
et le rayon de la trajectoire de S, puis en d^duire le 
rayon de la trajectoire de I. Cette solution est 
toutefois mo ins robuste que celle exposee ici, car le 
point S etant plus ^loign6 que le point I de I'axe de 
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rotation, les informations sur la position de S sont 
plus entachees d'erreur que celles sur la position de I. 
On utilise done les informations "siires" (celles sur I) 
comme base de la m^thode et non pas celles sur S, qui 
5 sont plus incertaines* Une autre solution consisterait 
egalement a trouver de mani^re simultanSe les rayons et 
centres des trajectoires de I et de S, Pour les m§mes 
raisons de robustesse, cette methode n'a pas ^t6 
retenue . 

10 Comment trouver les coordonn§es de I'isocentre, 

c'est-^-dire du centre de la trajectoire de I? On 
d^finit un plan tt comme etant orthogonal ^ Z' et ayant 
une cote suivant Z' egale d la moyenne des cdtes des 
points In* On projette orthogonalement sxir ce plan, en 

15 des points In tous les points . La trajectoire de I 

^tant suppos6e circulaire, les I»n sont a priori situes 
star un cercle* On cherche alors le cercle le plus 
proche, au sens des moindres carres, de 1' ensemble des 
I'n nous fournit: 

20 - le rayon ri de la trajectoire de I (c'est le rayon 

du cercle que 1 • on a trouv€) , 

- les coordonnees de I'isocentre, qui est ici defini 
comme 1 " intersection de n et de I'axe de rotation 
(c'est le centre du cercle que l*on a trouve) • 
25 Comment definir et trouver le cercle le plus proche? 
Soit C(A,r) un cercle de centre A et de rayon r. Ce 
qu'il faut, c'est definir une distance d(A,r) entre le 
cercle C(A,r) et 1' ensemble des I*n- Plusieurs m^thodes 
sont possibles. 
30 Methode 1 - figure 5a 

Soient (x'nfY'n) coordonnees de I'n- 

definit: d(A,r) = 2 I'n^'n 1^ point J'n etant defini 
comme 1 • intersection du cercle C(A,r) et de la droite 
Al'n- 
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Mdthode 2 - figure 5b 

Quand la mactiine de reconstruction 3D emet des tirs 
X suivant des positions angulaires r€gulieres, c'est 
dire tous les ^9 degree (SO constant) , on peut utiliser 
5 cette connaissance pour etre plus robuste dans la 
determination des paramdtres du cercle* 

on va dMinir d(A,r) par: d(A,r) = S I«n J'n/ 1^ 
point J'n etEint d^fini cette fois-ci comae le point du 
cercle C(A,r) qui possede, par rapport a une direction 
10 de reference, un angle polaire: e = Gq + (n-1) * SQ. 

Que l"on adopte la m^ttiode 1 ou la methode 2 pour 
def inir la distance entre le cercle et le nuage de 
points, la suite est la meme : on recherclie, psur une 
metliode quelconque de minimisation (gradient conjugu^ 
15 par exemple) les parametres A et r du cercle (ou,- pour 
la methode 2, les paramdtres a, r et 60) qui minimisent 
la distance au cercle* Une fois trouve ces pareundtres, 
c'est fini, puisque : 

- Le rayon de la trajectoire de J est connu : c'est r 
20 - L'axe Z» de rotation est entidrement connu, puisque 
I'on a d6termin6 son orientation a I'etape precedent et 
que l*on salt maintenant que cet axe passe par A. 
3* Trouver le rayon de la trajectoire de S 
On connalt le centre de la trajectoire de S, 
25 puisque c'est par definition le meme que celui de la 
trajectoire de I, que I'on vient de trouver a l»€tape 
pr^c^dente. Cherchons maintenant le rayon. On considdre 
les points S»n/ projections des points S sur le plem tt* 
Le rayon de la trajectoire de S est d^fini comme la 
30 moyenne des distances entre I'isocentre O et les points • 
S'n : 

rayon de I = l/N* SOS'n 

4* Trouver la difference de cate suivant l<axe Z» 
entre la trajectoire de la source et celle de I. 
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C'est la moyenne des cotes suivant I'axe Z' des 
points Sij: voir figure 5c. On en deduit la difference de 
cote par la distance au plan tt. 
5. Trouver 1» angle e: 
5 On travaille a partir des projections I' et S»n 

dans le plan tt. Pour chaque n, on calcule les angles 
polaires (par rapport a O) de I'n S'n on appelle 
en leur difference (voir figure 5d) . On estime e comme 
etant la moyenne des 6^* 

10 6. Trouver les tilts et 0£2 du d^tecteur 

Pour chaque vue, on calcule dans le rep^re 116 a 
l"axe de I'h^lice, les coordonn^es de la normale au 
detecteur. Ensuite, par application d'une formule de 
chcLngement de rep^re, entre le repdre lie a I'cixe de 

15 I'h^lice et celui li€ a I'isocentre et a l»axe de 
rotation Z* , on en deduit les coordonnees de la normale 
au detectexir dans le repere lie a I'isocentre et a l*axe 
Z" de rotation. Cette formule de changement de repdre 
est d^duite de la connaissance du repere lie a I'h^lice. 

20 Par ailleurs on connalt I'isocentre puisque c'est le 
centre du cercle trouve prec^demment . On connalt 
egalement 1 ' orientation de 1 » axe Z ' . Les deux autres 
axes de ce repere sont figes a une rotation prds. De 
preference on fait en sorte, pour les figer, que la 

25 premiere vue soit affectee dans le nouveau repdre d'un 
angle eo li^ a une position privilegiee de la machine. 
Dans ce nouveau repere, ces coordonnees ne dependent que 
des 2 angles et Pour chaque n, on calcule ces 2 

angles notes axxi a2n- Ensuite ai est d^fini comme la 

3 0 moyenne des am ®t idem poxir a2n* 
7. Trouver 1' angle at3 

Pour chaque vue, on calcule dans le repere lie a 
I'axe de I'helice, les coordonnees du vecteur unitaire 
parallele aux colonnes du detecteur. Ensuite, par 
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application de la formule de changement de rep^re, on en 
deduit ses coordonnees dans le repere lie a l*isocentre 
et a l"axe Z' de rotation. Or, dans ce nouveau repdre, 
ces coordonnees ne dependent que de 1' angle ^3. Poxir 
5 cliaq[ue n on calcule cet angle, note Q£3n^ Ensuite 03 est 
defini coiome la moyenne des a3n- 
Troisifeme etape 

On connalt maintenant 12 parametres intrinseques 
qui definissent la geometrie du systdme. A partir de la 

10 connaissance de ces 12 parametres intrinseques, il est 
tres facile de recalculer, pour chaque vue n, les 9 
parametres locaux de calibration def inissant cette vue. 
Ce calcul est d'autant plus simple qu'il a lieu dans le 
repere intrinseque de la machine, c'est-a-dire dans le 

15 repere lie a I'isocentre et a I'axe de rotation. On 
appelle (rsn, asn, zs^, thetan, phin/ psin, rin, ain, 
ziix) les 9 parametres locaux relatifs a la vue n, tels 
qu'ils sont definis dans 1» article publie au Congres de 
Rennes. II vient alors : 

20 rsn = rayon de la trajectoire de la source 

(= constante) 

asn = ©0 + iSO + (n-1) * Se 

zsn = difference de cote calculee precedemment 
(= constante) 
25 thetan = 9 0"*- (n-1) * 59 + al 

phin = a2 
psin = a3 

rin = rayon de la trajectoire de S (= constante) 
ain =90+ (n-1) * SQ 
3 0 zin = 0 (= constante) 

Ces nexif parametres ou coefficients de calibration 
sont done deteznaines a partir de 1' incidence deduite de 
9o + (n— 1) SQ» lis sont connus avec beaucoup de 
precision compte tenu des evaluations statist iques dont 



26 

ils r^sultent. Avec des mesures precises de 1» incidence 
on arrive alors a des coefficients de calibration 
calcules tres precis. L' incidence elle meme est deduite 
de I'indice n* Cet indice est connu par la division de 
la circonf erence de la machine en un nombre donne de 
positions pr€€tcLblies • 

On notera que, de preference, on recalculera, H 
partir des expressions analytiques selon les 
coefficients intrinsdques, les coefficients de 
calibration relatifs aux N vues calibrees, De ce fait, 
les coefficients de calibration pour ces N (100) vues 
sont ^valu^s avec autant de precision que po\ir les 
autres vues (les 900 autres) . Le resultat d* ensemble est 
encore meillexir. 

Une autre m^thode possible consisterait estimer 
directement les psirametres globavix de la machine 
partir de 1* ensemble des coordonn6es des billes 
d^tect^es dans toutes les projections de I'helice. Mais 
cette m€thode est longue. 
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REVENDICATIONS 

1 - Systdme d* imager le par rayons X avec un 
detecteur bidimensionnel sitiue en vis H vis d'un t:\ibe a 
rayons X, ce sys-beme d* imager ie par rayons X comportant 
un dispositif de calibration geom€trique miini d'un 

5 fantome de dimensions connues et de moyens de mesure 
pour d^duire, d pcortir d' images en projection du fantdme 
sur le d^tecteiir bidimensionnel, des coefficients de 
calibration de ce systeme d"imagerie par rayons X, 
cciract^rise en ce c[ue le fantdme comporte une 

10 succession, de proclie en proche, de structiires 
cell\ilaires a absorption radiologique contrast^e peur 
rapport d leur environnement , chague structure 
cellulaire §tant identifiable automatiguement par une 
caract^ristigue ordonn€e d • identification de cette 

15 structure cellulaire. 

2 - Systdme d'imagerie par rayons X selon la 
revendication 1, caract^rise en ce gue, ce detecteiir et 
ce tube toumant ensemble autour d'un axe, le fantdme 
comporte ^galement un axe, un support, et une 

20 distribution des structures cellulaires sur ce support, 
de telle fa9on gue des altitudes de ces structures 
cellulaires, mesur^es selon cet axe du fantdme, soient 
monotones et croissantes avec un num^ro d'ordre de ces 
structxires cellulaires, l»axe du fantome 6tant destine a 

25 atre orients dcins le syst^e d'imagerie par rayons X, en 
phase de calibration, sensiblement parall^lement ^ l"axe 
de ce systeme d'imagerie par rayons X, et de telle fa9on 
gue cette relation d'ordre se retrouve dans la 
projection du fantdme. 

30 3 - Systeme d'imagerie par rayons x selon la 
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revendication 2, caxacterise en ce c[ue le support est vin 
cylindre, de pr§f§rence circulaire, dont les 
generatrices sont pcuralldles ^ I'axe de ce cylindre et 
en ce que la succession ordonn€e des structxires 
5 cellulaires sur ce support a sensiblement la forme d'une 
he lice. 

4 - Systeme d'imagerie par rayons X selon la 
revendication 3, caract^rise en ce que le support 
cylindrique est muni a ses extr^mites de plaques de 

10 calage po\ir empScher sa deformation. 

5 - Syst&ae d'imagerie par rayons X selon I'une des 
revendications 1 ^ 4, caracteris^ en ce que la 
caract^ristique ordonn^e d • identification de chaque 
structure cellulaire est constitute par un reptrage de 

15 cette cellule par rapport au moins une structure 
cellulaire particuli^re du fantdme. 

6 - Systeme d'imagerie par rayons X selon la 
revendication 5, caract^ist en ce qu^une structure 
cellulaire particuli^re du fantdme possede une-taille ou 

2 0 xane absorption radiologique difftrente de celle des 
autres structures cellulaires et est placee sensiblement 
au milieu de la distribution des structures cellulaires* 

7 - Systeme d'imagerie par rayons X selon I'une des 
revendications l a 6, caracterist en ce que les 

25 structxires cellulaires possedent des absorptions 
radiologiques insuf f isantes pour empecber totalement le 
passage du rayonnement X a travers elles. 

8 - Syst^e d» imager ie par rayons X selon I'une des 
revendications 1 a 7, caracterise en ce que les 

30 structures cellulaires comportent des billes mttalliques 
plactes au fond de puits r^lises dans xin support. 

9 - Precede de calibration gtometrique d'un systtoe 
d"imagerie par rayons X a detecteur bidimensionnel 
tournamt autour d'xm axe, en vis ^ vis d'un tube ^ rayon 
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X et cozaportant les e1:apes suivantes: 

- on place un fan'tome de dimensions connues entre le 
tiibe et le d^tecteur, 

- on mesiire, pour une position en rotation de 
5 1" ensemble tube d€tecte\ir par rapport au fantome, dans 

1' image radiologique du fantome projet^e sur le 
detect e\ir bidimensionnel, des coordonn6es de lieux 
images de points caracteristiques du fantome, 

- on en deduit, pour cette position en rotation de 
10 1» ensemble txibe d^tecteur, et dans un rep^re associ^ a 

la position du fantdme, des coefficients de 
calibration relatifs aux positions respectives d'un 
foyer de rayonnement du tube a rayons X et du 
detectexzr bidimensionnel, 
15 - et on r^itdre ces deux derni^res Stapes pour des 

positions dSsirees differentes de 1' ensemble tube 
d€tecteur, 
caractSrise en ce que 

- on traite les coefficients de calibration relatifs a 
20 un certain nombre de positions en rotation de 

1* ensemble tube d€tecteur pour en extraire des 
param^tres intrinsdques (a^, a2, as . • • ) du systdme 
d'imagerie par rayons X qui sont independants de la 
position en rotation de 1' ensemble tube dStecteiir, 
25 - et on deduit de ces parametres intrinsdques des 

coefficients de calibration (rs^/ asn ) relatifs a 

une position quelconque en rotation de 1' ensemble tube 
detecteur en evaluant les valeurs de fonctions 
analyticjues de ces paramdtres intrinseques pour une 
30 valeur angulaire de position en rotation de 1" ensemble 

tube dStecte\ir. 

10 - Precede selon la revendication 9, caractSrisS 
en ce c[ue pour dedudLre des coefficients de calibration 
relatifs a une position desiree en rotation de 
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1' ensemble tube detectevur, 

- a) on simule, avec des coefficients de calibration 
arbitraires, des coordonn^es de lieux tb^oriques qui 
soient des images,, sur le detecteur bidimensionnel, de 
points caract^ristiques du fantdme, 

- b) on compcure les coordonnees de ces lieux 
theoriques avec celles des lieux mesurSs dans 1' image 
du fantdme, 

- c) on modifie en correspondance les coefficients de 
calibration arbitraires, 

- d) et on r^ittoe les deux stapes a) et b) jusc[u"d. ce 
que la comparaison montre un difference n^gligeable, 
auc[uel cas on retient comme coefficient de calibration 
poiar cette position en rotation de 1* ensemble tube 
detecteur les derniers coefficients de calibration 
modifies. 

11 - Proc^d^ selon la revendication 10, caracterise 
en ce que 

- on remplace les coefficients de calibration 
arbitraires par des coefficients de calibration 
deduits de la construction mecanique du systdme 
d* imager ie par rayons X. 

12 - Proced^ selon la revendication 10 ou la 
revendication 11, caract6ris6 en ce que 

- pour d^duire des coefficients de calibration 
relatifs d. une nouvelle position d6sir§e en rotation 
de 1* ensemble tube detecteur, d§cal^e par rapport k 
une position en rotation pr^c^dente, on utilise el la 
place des coefficients de calibration arbitraires 
relatifs a cette nouvelle position en rotation, les 
derniers coefficients modifies obtenus poiir cette 
position en rotation precedente. 

13 - Proc€d6 selon I'une des revendications 9 a 12, 
caract^ris^ en ce que les coefficients de calibration 
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sont: exprimes par des fonctions analytiques des 
paramdtres intrinseques, ces fonctions analytiques §tant 
d€crites dans iin repere cylindrique. 

14 - Proc6d€ selon I'une quelconque des 
revendications 9 ^ 13 caracterise en ce que les 
paxrametres intrinsdques sent obtenus par Evaluation 
statistiques a partir d'\an nombre restreint d' images du 
fantdme pour diverses orientation du systeme d« imager ie, 

15 - Procede selon l'\ine quelconque des 
revendications 9 d 14 caracterise en ce que on dSduit, ^ 
partir des parametres intrins^ques^ des coefficients de 
calibration relatifs ^ des positions en rotation de 
1' ensemble tube detecteur pour lesquelles on a effectue 
des mesures des coordoim^es de lieux images de points 
caracteristiques du fantdme, et pour lesquels on en 
avait d6duit des coefficients de calibration. 
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